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Abstract We consider a Brownian motion (BM) x { t ) and its maximal value Xmax = niaxo<T<t x { t ) 
on a fixed time interval [0, t]. We study functionals of the maximum of the BM, of the form Omaxit) = 
Jq V (xniax ~ x{T))dT where V(x) can be any arbitrary function and develop various analytical tools to 
compute their statistical properties. These tools rely in particular on (i) a “counting paths” method and 
(ii) a path-integral approach. In particular, we focus on the case where V{x) = S{x — r), with r a real 
parameter, which is relevant to study the density of near-extreme values of the BM (the so called density 
of states), p{r,t), which is the local time of the BM spent at given distance r from the maximum. We 
also provide a thorough analysis of the family of functionals Ta{t) = (a^max—a;(r))“ dr, corresponding 
to V(x) = x°‘, with a real. As a is varied, Ta{t) interpolates between different interesting observables. 
For instance, for a = 1, Ta=i{t) is a random variable of the “area”, or “Airy”, type while for a = —1/2 it 
corresponds to the maximum time spent by a ballistic particle through a Brownian random potential. 
On the other hand, for a = —1, it corresponds to the cost of the optimal algorithm to find the maximum 
of a discrete random walk, proposed by Odlyzko. We revisit here, using tools of theoretical physics, 
the statistical properties of this algorithm which had been studied before using probabilistic methods. 
Finally, we extend our methods to constrained BM, including in particular the Brownian bridge, i.e., 
the Brownian motion starting and ending at the origin. 
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1 Introduction 

Stochastic processes are at the heart of many fundamental problems in statistical physics. In particular, 
it was realized a long time ago that Brownian motion (BM) is the process underlying many physical 
systems and corresponding random models. Since then, BM has not only become a cornerstone of 
statistical physics [TIEIIH] but has also found numerous applications in various areas of science, including 
biology j4], computer science [5l|6j[7l[8l|9] or financial mathematics mii\- These various applications 
have motivated the study of functionals of Brownian motion |Hl lTni[T^ll8| . which are observables of 
the form 0{t) = V(x(t)) dr where V(x) can be any function and x(t) is a BM (see Fig. 1 (a)). 

It might also be relevant to consider functionals of variants of BM, like the Brownian bridge (BB), 
xb{t), which is a BM conditioned to start and end at the origin (see Fig. (b)), the Brownian 
excursion (BE), xe{t), which is a BB conditioned to stay positive on the whole time interval [0,t] 
(see Fig. [b (c)) as well as the Brownian meander, XMe('r), which is constrained to stay positive on 
[0,t] but can end up at any point at time t (see Fig. [^(d)). For instance, if V{x) = S(x — a), the 
Brownian functional 0{t) corresponds to the local time at the fixed level a, which is an important 
quantity in probability theory |14| . Another interesting example concerns the case where V(x) = x for 
a BE, which corresponds to the area under a BE. In this case, the distribution of Osit) = /n ^ e{t) dr 
is given by the so called Airy-distribution which appears in computer science |15l[THlfr7l[Ib| as well 
as in the extreme statistics of elastic interfaces EiiaiiD]. Extensions of the Airy-distribution to the 
area under Bessel processes {i.e., radius of the d-dimensional process) have been recently discussed in 
Ref. |21| . Yet another example which is relevant in finance is the case where V{x) = exp (cc), which 
describes the price of an Asian stock option in the Black-Scholes framework |^. Note that in this 
exponential case, 0{t) = exp (x(t)) dr also represents the stationary current of a disordered Sinai 
chain connected to two reservoirs of particles Quite interestingly, these functionals of 

BM and its variants can be studied using powerful tools of theoretical physics, namely path integrals 
methods (leading to the so called Feynman-Kac formula). This allows to recast the study of Brownian 
functionals in a quantum mechanical framework Pll?7| . 
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Fig. 1 Brownian motion x(t) (a) and its variants: the Brownian bridge xb{t) in (b), the Brownian excursion 
xe{t) in (c) and the Brownian meander XMeir) in (d). 


While such Brownian functionals are thus well understood, several recent works have pointed 
out the physical relevance of functionals of the maximum of Brownian motion, which are much less 
studied. In this case, one is interested in observables of the form Omax(t) = /q M(a;max — a;(r))dr 
with a;niax = inaxo<T<t a^(T) where x{t) is BM or one of its variants (see Fig. [R. An important case 
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corresponds to the case where V{x) = S{x — r) where 

p{r,t)= / i5(xi„ax - a;(T) - r)dr , ( 1 ) 

Jo 

which is the so called density of states (DOS) near the maximum. This is a natural and useful quantity 
to characterize the crowding of near-extremes |28| . Indeed, p{r^t)dr denotes the amount of time spent 
by x{t) at a distance within the interval [r,r -\- dr] from x^aa^ (see Fig. |^. Hence p{r,t) is similar to 
the local time with the major difference that here the distances are measured from ccmaxi which is itself 
a random variable. The statistics of the DOS was recently studied by us in the context of near-extreme 
statistics |1S]. Note that, by definition, p{r,t) dr = t. Therefore its average value, {p{r,t))/t, where 
(...) means an average over the trajectories of BM has a natural probabilistic interpretation as it is 
the probability density function to find the BM at a given distance r from the maximum in the time 
interval [0,t]. In particular, the average value of any functional of the maximum can be expressed as 

(Omax(O) = ( / V{Xma^ - x{t)) dr) = / {p{r,t))V(v) dr , ( 2 ) 

Jo Jo 

which naturally holds not only for BM but also for its variants, like Brownian bridge or Brownian 
excursion. 



Fig. 2 (Color online) One realization of the process x{t) on the time interval [0,t], with a width W{t) = 
max.rg[o,ti x{t) — min,-g[o_t] x(t). x(r) spends a time p(r, t)dr at a distance within [r, r -\- dr] (the green stripe) 
from the maximum Xmax, with p{r,t) being the DOS 0 - 


Another very interesting application of functionals of the maximum of Brownian motion concerns 
the case where V{x) = l/(2a:), which enters into the analysis of the optimal algorithm to find the 
maximum of a discrete random walk of n steps with n ^ 1 m- Indeed, let us consider a discrete 
random walk (RW), starting from Aq = 0 and evolving via the Markov rule: Xk = X^-i + rjk where 
77 fc = ±1 with equal probability 1/2. We study the search problem of finding the maximum of the RW, 
M„ = maxo<i<n Ai, while minimizing the number of values of A^’s that are probed. The cost of the 
algorithm is identified by the number of probes used to find Mn- The simplest algorithm consists of 
probing all positions A^ for 1 < fc < n: its cost is n. In Fig. we show an example for n = 14 where 
we find M 14 = 7 in 4 probes: this shows that, because of strong correlations between the positions of 
the RW, n is actually a rough upper bound of the cost of the optimal algorithm. Of course, some RWs 
need more probes than others in order to find M„. For example, the RW with n jumps 4-1 needs only 
one probe in A„ = n but the RW with alternating jumps ±1 {X 2 k = 0 and A 2 fe+i = 1) needs -|- 1 
probes to be sure that none X 2 k = 2. But of course these cases are rare. Let An be the ensemble of 
algorithms that find Mn, according to the above rules, and let us denote by C{a) the cost, as defined 
above, of the algorithm a that belongs to An- Of course C{a) is a random variable, which varies from 
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one realization of the RW to another, its average value being denoted by {C{a)). In Ref. [3^, Odlyzko 
studied the minimal average cost of such algorithms and he showed that, for large n, the minimal 
average cost is proportional to i/ri, much smaller than the cost of the aforementioned naive algorithm, 
which necessitates n probes. Indeed, one has EHl 

min {C(a)) = Co^/n + o(y/n) , (3) 

a^A„ 


where cq is a constant given by 


CO = (/),/ 


1 dr 

2 Jo [Xmax - x{t)] 


(4) 


hence the relevance of the functional of the maximum with V{x) = l/{2x) for this search problem. 
Note that Odlyzko found an expression of cq in terms of a complicated double integral which was then 
evaluated independently by Hwang m and Chassaing [32] to be 


Co 



log2 = 1.1061... . 


(5) 


Furthermore, Odlyzko found an algorithm, denoted as Od(n), belonging to An, such that 


(C'(Od(n))) = coV^ + o(i/n) . 


( 6 ) 


Hence Eqs. (§ and (§ show that Odlyzko’s algorithm Od(n) is on average asymptotically optimal. 
In a subsequent work, Chassaing, Marckert and Yor showed that Od(n) is not only asymptotically 
optimal on average (|^ but also in distribution [33], which means the following. If one defines d>n{x) 
as 

d>n{x) = min Pr (, (7) 

y y'n J 

then one has for any x |33| : 


lim <Pn{x) = Pr(/ > x) , (8) 


where I is the random variable defined above in In Ref. [33], the authors studied the distribution 
of / as well as its moments, using rather involved probabilistic methods. We will show here how these 
results can be derived simply using path integrals techniques. 

Yet another case of a functional of the maximum, Oynax{t), corresponds to the case where V{x) oc 
l/v^. In this case Omax(0 describes the largest exit time of a particle, of unit mass, moving ballistically 
through a random potential on the segment [0,1] (at zero temperature). Consider indeed a random 
Brownian potential x{y) over a line segment y S [0,1]. Imagine shooting a classical particle of unit 
mass with fixed energy E from the left of the segment at y = 0. The energy conservation leads to 


1 

2 



+ x{y) = E . 


(9) 


Clearly, this classical particle can penetrate the region y G [0,1] if and only if its energy E is bigger than 
the maximum value of the potential x{y) over y G [0,1], i.e., HE > ccmax where Xmax = niaxo<j,<i [a;(y)]. 
Now, imagine sending a beam of classical particles with varying energy through this potential barrier. 
Only those particles with energy larger than Xmax will go through the barrier. The time taken for such 
a penetrating particle to exit the region y G [0,1] through its right can be computed from Eq. ([^ as 


T{E) 


J_ dy 

Jo \/E - x(y) ’ 


( 10 ) 


where E > Xmax- Clearly, the time needed by a penetrating particle to cross the region y G [0,1] 
depends on the energy E of the particle. The slowest particle, i.e., the one that takes the longest time 
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to cross, is the one that has the lowest allowed energy to penetrate, i.e., the one with energy E = x„ 
Hence, the maximum time needed by a particle to cross the barrier is given by 


^Tmax — 


1 

inaxE>x^,AT{E)] = [T(a;max)] = ~^ 


dy 


\/ a^max - x{y) 


( 11 ) 


which thus corresponds to a functional of the maximum of BM, Omax(i = 1), with V{x) = In 

view of these two interesting physical examples in Eqs. and ( |11[ ), it is rather natural to consider 
the family of functionals of the maximum of BM Oj^axit) with V(^ = a;“ such that 


Ta{t) = / dr (a;niax - x{t)Y 

Jo 


( 12 ) 


indexed by a real a €] — 2, oo[. Note that, by using the self-affinity of BM, Ta{t) = 1). 

In particular for a = — 1 this corresponds to I in 0 while for a = —1/2 this corresponds to Tmax in 
([Tl] ). On the other hand, for a = 1, is called the area under a Brownian double meander |18l 

120] . For other values of a, Tait) generalizes these three cases. 

Of course, one can consider similar observables as in (121 for the Brownian bridge, i.e. when 
x{t) —>■ xb(t) (see Fig. [^(b)) and Xmax.B denotes its maximum on the time interval [0,t]: 




Jo 


Xmax.B “ XB(r))“dr 


(13) 


with, as above, Ta^sit) = !)• The simplest case is a = 1, which corresponds to the 

area under a Brownian excursion: this can be easily seen by permuting the pre-minimum and the 
post-minimum part of a Brownian bridge, i.e. by using Vervaat’s construction |36| . Hence ,B{t) 
is distributed according to the Airy distribution, discussed above 1191120] . For a = —1/2, this yields 
again the maximal time spent by particles to pass through a disordered periodic potential which is 
a Brownian bridge. Finally we conjecture, following the lines of reasoning of Refs. [soiisa, that for 
a = —1, i. e. the equivalent of / in (ffl where x{t) is replaced by xb{t) and a^niax by Xmax.Si ^(t) 

yields the cost of the optimal algorithm to find the maximum of a random walk bridge (see Appendix 


C.3). As shown below, this random variable turns out to be related to the maximum of a Brownian 


excursion (see also Refs. 131135]). 


2 Summary of main results 

The goal of this work is to present various tools to study the statistics of such functionals (!lmax(0 of 
the maximum of Brownian motion and its variants. It is useful to summarize the different approaches 
developed here as well as the main results obtained in the present paper. The first natural observable 
to compute is the average value of such functionals, (Omax(f)). According to Eq. the average 
can be obtained by computing the average density of states (p(r, <)). As we show here, this quantity 
can be calculated using rather elementary computations, based on the propagator of the Brownian 
motion with appropriate boundary conditions. In particular, this method allows us to recover in a 
very simple way the result for cq given in Eq. 0 and previously obtained in Refs. Ell and [32] from 
the analysis of a rather complicated double integral. The expression (p(r, t)) for the free BM, as well 
as for the bridge, was recently announced by us in a short Letter |29] . In this paper, in addition to 
providing details of these computations, we extend our techniques in several directions obtaining many 
new results. In particular, we show that this method, relying on propagators, can be easily adapted to 
a variety of other constrained Brownian motions, including the excursion, the meander, as well as the 
reflected BM and the reflected Brownian bridge. The main characteristics of {p{r,t)) for these various 
constrained BM are summarized in Tables and below. From it, we obtain in particular from (|^ 
the average of the functional Ta{t) for the free BM, with a g] — 2, -|-oo[ as 

(2t)i+t(2-2-“)T(i^) 

(2 -I- a)^/7T 


(Ta(t)) 


(14) 
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from which we get in particular the constant Cq in Eq. (j^ as cg = Ta=-i(t = l)/2 = -y/ 8 / 7 rlog 2 = 
1.1061... as given in Eq. pi. For a = 1, one recovers = 1)) = [50] for the average 

area under a Brownian double meander (see Fig. [^. As a function of a it has an interesting non¬ 
monotonic behavior, diverging when a —>■ —2 as {Ta{t = 1)) ~ 4/(2 -|- a) as well as when a —)■ oo as 
{Ta(t = 1)) ~ ja, exhibiting a minimum for a ~ 1.148. Similarly, for the bridge one 

obtains 


A +“/2 / Q,\ 

= 2) ■ (15) 

In particular, for a = 1 it yields back the first moment of the Airy-distribution, {Ta_=i ,B{t = 1 )) = 
y^Tr/S [371138] while for a = —1 it gives the equivalent of cg in Eq. (|^ for the Bridge, cq^b = 
Ta=-i,B(t = l)/2 = •\/ 7 r /2 = 1.25331.... Interestingly, this means that, on average, the cost of 
Odlyzko’s algorithm is higher for the BB than for the free BM. This can be roughly understood 
through the fact that the DOS close to the maximum is slightly higher for the BB, which is pinned to 
the origin on both sides of the time interval, than for the BM which is free on one side. As a function 
of a, (T„, B{t)) is also non-monotonic diverging when a —>■ —2 as {T^^Bit = 1)) ^ 4/(2 -|- a) as well as 
when a —>■ 00 as {Ta^Bit — 1)) i/tto, with a minimum for a « 2.960. 


Brownian Motion (BM) 



a 


Brownian Bridge (BB) 



a 


Fig. 3 Left panel (bottom): Plo t of the variance Var(ra(t = 1 )) = {T^{t = 1 )) — {Ta{t = 1 ))^ as a function 
of a for the BM, as given in Eq. ( | 16 [ ). Left panel (top): Focus on the asymptotic behaviors of Var(ra(t = I)) 
for a —>■ —2 and a —>■ 00 . The dasUed lines indicate the asymptotic behaviors given in the text below Eq. (161. 
Right panel (bottom): same quantity for the Brownian bridge, Var(Tc,s(t = 1 )), as given in Eq. (T7l. 
Right panel (top): Focus on the asymptotic behaviors of Var(rci,s(t = 11 ) fo r a —^ —2 and a —>■ 00 . the 
dashed lines indicate the asymptotic behaviors given in the text below Eq. ( | 17 [ ). The functionals Ta{t) and 
Ta,B{t) interpolate between different interesting observables: for a = —1 they describe the cost of Odlyzko’s 
optimal algorithm to find the maximum of a discrete random walk, for a = —1/2, they describe the largest exit 
time of a particle through a Brownian potential and for a = 1, Ta=i{t) and Ta=i,B{t) are random variables 
of the “area” type. For the bridge (right panel), this corresponds to the area under the BE (Airy random 
variable). 


If one is interested not only in the average of the functional Ta{t) but also in higher cumulants 
or even in its full distribution, the knowledge of the average DOS is obviously not enough and these 
computations based on propagators become cumbersome. Instead, we present here a method, based 
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on path-integral which allows us to compute its Laplace transform, namely from which 

the moments of arbitrary order can be obtained by differentiation with respect to (wrt) A. In some 
cases, this can also allow us to compute the full distribution of the functional. We first illustrate 
this method on the DOS itself, for which V{x) = 5{x — r). The main results in this case were 
recently announced in Ref. |29| . without any details, which we thus provide here. Then we use this 
general formalism to study the special functionals Ta{t) (12 1 as well as Ta^sit) in (ll^ for the BB. In 


particular, we show how this method allows us to compute the second moment, and eventually the 
variance YaiiTait)) = — {Ta(t))'^ for arbitrary a g] — 2 , -|-oo[ under the form: 


Var(T,(t)) = 


^ 2 +a 


(2“ - l)(2“+i - l)T(a -h 3) 22“+2(2 “+i - 


(23“(q;-I- 2)2 y (a-f 1)2 tt 

{a + 2)r{2a + 2)(22“+2 - 2 “+i(a -f l)Ri/2(a + 2 , - 2 (a -f 1 )) - 1 ) 
^ (a + l) 2 T(a + l) 


(16) 


where ^^(a, h) = fj ^(1 —t)*” ^ dt is the incomplete beta function. Formula (16) yields in particular 

Vai(Ta=-i(t = 1)) = ^ + (4 — 32/7r) log(2)^, as obtained previously in Ref. [33], using probabilistic 
methods. For 0 = 1, Eq. yields Var(Ta-i(t = 1)) = 17/24 — 2 / 7 r |20]. Interestingly, as a function 
of a, it has a non monotonic behavior. It is diverging when a —?■ —2 as Yai(Ta(t = 1)) ~ 8/(q; -I- 2) as 
well as when a —oo as 4\/2{^a je)'^o~'^ exhibiting a single minimum for a = 0. In Fig. [^(left panel), 
we show a plot of Var(To,(t = 1)) as a function of a. 

Similarly, one can also compute the second moment in the case of the bridge, yielding the variance 
Var(r„,B(f)) = - (ra,B(t))2: 


Var(r„,B(t)) 


V^(T(2a-f 3) -T(a-h2)2) 
(a + l) 2 T(a+ |)23“+i 


n§ + l)^ ' 

2 “ 


(17) 


In particular, we can check that Yai{Ta=i^B{t = 1)) = 5/12 — tt/S which is the variance of the Airy 
distribution [m[3ini371l3H] . while for a = —i, one has Var(TQ,^_i ^(t = 1)) = 27r^/3 —27r. As a function 
of a it has also a non-monotonic behavior, diverging when a —>■ —2 as YaiiTa^sit = 1 )) 8 /(a -|- 2 ) 

and when a —>■ oo as y^Sir/a, exhibiting a minimum for a = 0. In Fig. (right panel), we show 
a plot of Var(T„ B(t = 1)) as a function of a. 

Finally, in the special case a = —1, which corresponds to the cost of the Odlyzko’s algorithm, we 
are able to compute exactly the moments of arbitrary order, both for the free BM, = 1)) and 

for the bridge b(^ = !))• obtains indeed. 


(TL-i(t=l))=A 


fc + l\ 2 t +2 


k 

Y. 

m—0 


C(m)C(fc — m) , C(^) = (1 — 2 ^ ^)C,{m) = 


Y 

n>l 


(- 1 ) 


n+1 


,(18) 


recovering, using a completely different method, the result of Chassaing, Marckert and Yor [53]. For 
k = 1 and k = 2, this formula (18) yields back the results for the first moment (14) and for the 
variance (16). 


For the bridge, we obtain the result, for any real k 




1+1 


-fc- 


^kr -- 


ai-k). 


(19) 


Of course for fc = 1 and k = 2, this formul a (|l^ yields back the aforementioned results for the first 
moment (15) and the variance ( [T7| . In Ref. |33| . the authors obtained the full probability distribution 
function (PDF) of Ta=-i{t = Ijm terms of a convolution of two theta-functions. Here, we obtain the 
full PDF pb{s) of Ta=-i^B{t = 1) in the case of the BB, as: 


!>oW = *PE( 1 - 




7T1=0 



8 ^ 2713/2 


771=1 


( 20 ) 

























where the two formulas are related to each other via the Poisson summation formula. Interestingly, in 
Eq. (201, one actually recognizes the PDF of the maximum of a Brownian excursion on the unit time 


interval, a;niax,£; = niaxo<r<i xe(t). One has indeed 


Pb{s) = ^Pr < 0 ■ 


( 21 ) 


As explained below, one can show that Ta=-i ,b( 1) /o^ dT/a;£)(r). Hence, Eqs. J20|) and 

manifestation of a non-trivial identity in law for the Brownian excursion xe(j) IME 


dr law „ , . 

-= 2 max xe{t) . 

xe{t) 0<t<1 


( 22 ) 


Therefore, our result for Ta=-i^Bit) in Ed- ( |20[ ) provides a simple derivation of this non-trivial identity 
(221, which was proved in Refs. |Mll^ using rather involved probabilistic tools. 

The paper is organized as follows. In section]^ we focus on a method based on "counting paths", 
using propagators of BM. We first illustrate this approach to compute the average DOS for BM (p(r, t)) 
in subsection |3.1|and then extend it to the Brownian bridge in sectionas well as to other constrained 
BMs in sections |3.3| and [3)^ and Appendix |B| In section [SAl we also present a comparison between our 
exact results with numerical simulations, the details of which are given in Appendix]^ In section]^ we 
develop a general path integral approach to compute functionals of the maximum of Brownian motion 
and Brownian bridge. Within this framework, which is presented in detail in section |4.1| we study 
more specifically th e full s tatistics of the DOS in section [4.2.1 and the family of functi onal T a(t) [see 


Eq. (@1 in section [4.2.2[ We also analyze more precisely the case a = — 1 in section [4.2.3[ which is 


relevant to study Odlyzko’s algorithm. Some technical details, including a description of the main 
ideas behind Odlyzko’s algorithm, have been left in Appendices. 


3 The method of propagators 

In this section, we present a rather simple method, based on "counting paths", to compute the average 
DOS {p{r,t)). In the cases of BB and BE (which turn out to be identical as we show below), this 
method allows us also to compute higher moments {p%{r, t)) with k an arbitrary integer. We illustrate 
the method in detail for the case of free BM and then apply it to various constrained BM: the Brownian 
bridge and Brownian excursion in section 3.2 and the Brownian meander in section 3.3 In Appendix 


[^we use the method, for completeness, for the reflected BM and for the reflected BB. 


3.1 The average DOS for free Brownian motion 

To compute (p(r, t)), we simply average Eq. 0 over the trajectories of the BM and write it as 


(p(r,t))= / ((5(a:max - a;(r) - r)) dr . 


(23) 


In Eq. (231 we recognize that the integrand ((5(a:max~a^(r)—r)) has a simple probabilistic interpretation. 
Indeed, one has 


{S{xnia^ - x(t) - r))dr = Pr [x{t) G [a:max - r - dr,: 


(24) 


The idea to compute the PDF {S{xmax — x(t) — r)) in (241 is to evaluate the "number" of paths that 
reach their maximum M at time tmaxi pass through M — r at time r, and end at < M at time t. 
The total number which we want to compute is then obtained by integrating over xe, M and tmax (see 
Fig. g. In each of the three time intervals delimited by r, tmax and t, the BM is constrained to stay 
below M. This number of paths can thus be computed from the propagator GM{o!.\P,t) of the BM, 
starting at time 0 at cc = a < M and arriving at time t at x = /3 < M and staying below M during 
the whole time interval [0,t]. It can be easily computed using, for instance, the method of images: 


GM(a|/3,t) = 


1 


. CSzisA 




y/2TTt \ 


(25) 
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For future purpose, it is useful to compute its Laplace transform (LT) wrt t, Gm{ci:\P, s). It is given 
by 

r°° 1 / \ 

GM(a|/3, s)= t)dt=^ L-VTs\ii-^\ _ ^26) 

Jo v2s V / 

When dividing the time interval [0, t] into three parts as in Fig. two cases may arise: r < or 



Fig. 4 (Color online) Illustration of the method of propagators. The BM starts at 0 and visits the point M — r 
at time r, passes through M — e aX time tmax and ends at at time t. 


T > tmax- One can show that these two configurations give rise to the same contributions to {p{r,t)): 
this can be seen by making a global shift x{t) —>■ x{t) — xf/2 (for each realization of BM) and using 
the time reversal symmetry. Note also that, when dealing with BM which is continuous both in space 
and time, one can not impose simultaneously x(traax) = M and x(t) < M right before and after tmax- 
To circumvent this difficulty, one imposes instead x(tniax) = M — e, and take the limit e — > 0 at the end 
of the calculation. Finally, {p{r,t)) is given by the ratio of paths which have the analyzed properties 
to the normalization constant Z{e) which count the same paths without the condition to pass through 
M — r at time r 

2 pOO pt pJVl pim&x 

{p{r,t)) = lim—— / dM / dt^ax / dxp / dr Gm(0|M - r, r)GM(M - r|M - e, ti„ax - r) 

e^oZ(e) Jo Jo J_^ Jo 

xGm{M - e\xF,t - ( 27 ) 


where we have used the Markov property of BM and where the factor of 2 comes from the two 
aforementioned equivalent situations corresponding r < tmax and r > tmax- In (271 the normalization 
constant Z{e) is given by 


pOO nt pM 

Z{s)= / dM / dfinax / da:FGM(0|M-e,tmax)GM(M-e|xF,t-tmax) ■ 
Jo Jo J — oo 


(28) 


The normalization is easily computed as Z{£) ^ 2e^ when £ —>■ 0. In (271, we recognize a convolution 
structure. Taking the LT wrt to t, we find 


dte "*(/o(r,t)) = 


( 29 ) 


rM 


lim , 
e-yO Z{e) 


dM / dxF Gm{0\M — r,s)GM{M — r\M — e,s)GM{M — e\xF, s) . 
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Using (261 and performing a small e expansion, one obtains straightforwardly 

^—\/2sr /*^r 


I dte =8 / dM ^ smh{'/^M)e 

lo LJq U2s 

pM 

sinh (v^r)e-^’' / 


dxFC 


— y/^{M — X 


^OO -U2sM 

dM- 




Performing the remaining integrals over xp and M we obtain 

roo 

/ dte"'*‘(p(r,t)) = 8 


g-%/2sr _ g-2%/2sr 


(2s)3/2 


By inverting the above LT (311 we finally obtain 


(p(r,t)) = v^p , p{r) = 8 - <P^^\2r)^ , <P^'^\r) = 


- _!_-erfc ( — 1 

2 \V2) 


(30) 


(31) 


(32) 


where ^*^2) belongs to a useful hierarchy of functions, <P^^\ as explained in Appendix [a| From the 
average value of the DOS, one can compute the average value of any functional of the maximum, 
according to (|^. In particular, for the special case V{x) = x°‘, one obtains from ([^ and 


nt nOO 

{Ta{t)) = {J [Xmax - a;(r)]“dT) = j (p(r,t))r“dr 


(2 + a)y/7r 


(33) 


as announced in the introduction in Eq. (14l. In Fig. (left panel), we show a plot of {Ta{t)) as a 
function of a. Specifying this formula (331 to the case a = —1, one obtains cq = Ta=-i{t = l)/2 = 


log 2 as announced in (51, recovering in a rather simple way the result of Refs. mm- 


Brownian Motion (BM) 


Brownian Bridge (BB) 



6 



to 


a 


a 


Fig. 5 Left panel: Plot of {Ta{t = 1)), as a function of a, for the B M, a s given in Eq. ( |33[ ). Right panel: 
Plot of {Ta,B{t = 1)), as a function of a, for the BB, as given in Eq. ( |44| . In both panels, the dashed black 
lines indicate the asymptotic behaviors discussed in the introduction. In particular, for a = — 1 for the BM 
(corresponding to Odlyzko’s algorithm), we recover {Ta=-i{t = 1)) = 4^2/7rlog2 |31II32| . For q = 1 for the 
BB (corresponding to the Airy distribution), we recover (rct=i,fl(t = 1)) = \/vr/8 |19II20| . 
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This method can in principle be adapted to compute higher moments of the DOS, {p^{r,t)), but 
such computations involve a rather cumbersome combinatorial analysis, for fc > 1. A more powerful 
approach amounts instead to compute the generating function of p{r, t) - as shown in section 
As we show now, this method can also be easily extended to compute the average DOS for various 
constrained BM. 


3.2 The average DOS and its higher moments for Brownian bridges and Brownian excursions 


For Brownian bridges, as well as Brownian excursions, the method based on propagators of constrained 
BM, allows us to compute arbitrary moments of the DOS, {p%{r, t)) . An expression for these moments 
can be written from the definition [see Eq. Q where x{t) is replaced by the BB xb{t)] as 


{p%{r,t)) = / ... (Y\dtiS{xrni,^^B - XB{ti) - r)) 

Jo Jo 


(34) 


such that in (341 the BB visits k times the point M — r at successive times ^1,0, ■ • • ,^fc where M is 
the value of the maximum on [0,t]. To compute this quantity (341 for the BB, it is useful to invoke 
the Vervaat’s construction to relate the DOS of the Brownian bridge to the local time of the Brownian 
excursion (see Fig. [^. 


PB Tioc,E PE 



Fig. 6 (Color online) Illustration of the identities in law given in Eq. (35 I in the text relating the DOS 
near the maximum for the BB to the local time, on the one hand, and to the DOS, on the other hand, for 
the BE. (a) A typical trajectory of a BB, xb{t). The process spends a time pB{r,t)dr in the green stripe 
[a^max.B — r — dr, Xmax.s “ r]. (b) The process —y(r) = ®max,s — xb{t) after the Vervaat’s transformation 
explained in the text is an excursion. The time spent in the green stripe is now given by T\oc,E{'c,t)Ar for this 
new process y{T). (c) The process z{t) = 2 ;s(t) — Xmin.s after the transformation explained in the text is also 
an excursion. The time spent in the green stripe is now given by pE (r, t)dr for this new process. 


The Vervaat’s construction works as follows jM]: we denote by tmax the time at which the BB 
reaches its maximum on [0,t] (see Fig.[^(a)) and break the time interval [0,t] into two parts, [0,tniax] 
and [tmaxjt] as illustrated in Fig.|^(a). Then we permute the two associated portions of the path, the 
continuity of the path being guaranteed by XB{t) = Xb(0) = 0 for the BB. We can further transform 
the path by considering 2/(t) = Xmax.B — xb(t) and take finally the origin of times at f„iax: this yields 
the configuration shown in Fig.|^(b). Furthermore, if we initially break the time interval into [0,tinin] 
and [tminjt], where tmin denotes the time at which the BB reaches its minimum on [0,t], we permute 
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the two associated portions of the path and finally take the origin of times at tniim we obtain the 
configuration shown in Fig. (c). 

These transformations show that p£;(r, t) for an excursion is identical in law to the local time in r 
for an excursion, Tioc,£;(r, t) and hence also identical in law to the DOS of the BB. Hence the Vervaat’s 
construction allows to us write the following identities: 

PB{r,t)^= PE{r,t)^= Tioc,E{r,t) = [ S{xE{r) - r) dr , (35) 

where we emphasize that Xe{t) is a Brownian excursion. Therefore, the moments of the DOS for the 
bridge or the excursion can be computed as 


{P%{r,t)) = {p%{T,t)) 



(lU5(xE{ti) 


r)) . 


(36) 


As done before for the free BM in Eq. ( [2^ , the expression in ( [M| ) can be computed using the propagator 
of the excursion, i.e. the propagator of a free BM starting at a > 0 at time t = 0 and 

reaching /3 > 0 at time t and staying positive over the whole interval [0,t]. This propagator G~^{a\f3,t) 
can also be computed simply by the method of images, yielding: 


G+(a|/3,t) = 


1 








•\/27rf \ 

together with its LT transform wrt t, which will be useful in the following: 




,-V2s|a-/3| _ -V^ifi+a) 


)■ 


(37) 


(38) 


Since we are studying here a Brownian excursion, we need, as before, to introduce a regulator such 
that xe{ 0) = XE^t) = e > 0 and take the limit e —>■ 0 of a suitably defined quantity, as done in Eq. 
(27l - since we can not impose simultaneously xe{0) = 0 and a;£;(0''') > 0. One has here: 


{PBir.t)) = (p|(r,t)) = lini 


k\ 


e-i-O Zsie) 


/ dti---dtk G+(£|r,ti)G+(r|r,t2 - G) ■ 

'0<tl<t2< '<tfe<t 


■G+(r|£,t-4) , (39) 


where the combinatorial factor kl comes from the different permutations of the intermediate times 
ti, - ■ ■ ,tk where the Brownian excursion xe(j) reaches the value r. The denominator Zb{e) is given 
by 


ZB{e) = G+ie\e,t) 



£ —>■ 0 . 


(40) 


To compute the multiple integral over the times ,tk in the right hand side of Eq. 

recognize, as before for the free BM, a convolution structure. Hence its LT wrt t is given by: 


(391, we 


die dti---dtk G+(£|r, ti)G+(r|r, t2 - G) • ■ ■ G+(r|£, t - 4) = G+(£|r, s)^G+(r|r, s)'" ^ 

J0<ti<t2< -<tk<t 


, _ — „-2%/2sr'|/c-l 

Ap2 -2V^r\E_Z _ > 

(2s)G-i)/2 


£ —>■ 0 , 


(41) 


where we have used the explicit expression of the LT of the propagator in (38). Using the above 
expression (41) together with (40), one obtains finally after Laplace inversion: 


{p%{r,t = 1)) = (p|(r,t = 1)) = (T£„£;(r,t = 1)) = 2V^fc! ^(-l)'+i (2r0, (42) 




with the convention = 

(1421. Thus we recover in 


and where the are defined in Appendix [a| see Eq. 

the result obtained by Takacs [51] (see his Eq. (38)) by a probabilistic 
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method (note the correspondence between the functions J^’s in [33] and the functions Jfe(r) = 

/i°° da; ^‘^{x — 1)^ = 

For k = 1, one finds the mean DOS for the BB on the unit time interval, {psir^t = 1)) = 4re“^’’^, 
as found in Ref. [39]. Note that it coincides in this case with the PDF of the maximum of a BB, which 
is a generic property for periodic signals such that x{t) = a;(0) [ID]. One has indeed 


{pB{r,t = 1)) = {pE{r,t= 1)) 



dT(^(Xniax,i? 


— ((5[Xmax,i? ^]) 


XB{T)-r]) 


(43) 


where we have used the periodicity of the bridge and the possibility of adding to it an arbitrary 
constant. From the average DOS, we can compute, using (i), the average value of any functional of 
the maximum of a BB. In particular, for the interesting family of functionals V{x) = x°‘, one obtains 


pt pOO 

{Ta,B{t)) = { / (a^inax.B “ dr) = / {pB{r, t))r°‘ dr 

Jo Jo 


= t 


1+^ 


nl + f) 

2“/2 


(44) 


(45) 


as announced in the introduction in Eq. (15l. In Fig. (right panel), we show a plot of (Tq, ^(t)) as 
a function of a. 


3.3 The average DOS for the Brownian meander 


Using the same method based on propagators (271, we can also compute the average DOS for the 
Brownian meander (see Fig. (d)). To this purpose, we need to know the propagator of a Brownian 
particle confined in a given interval [0, M] with absorbing boundary conditions both in x = 0 and 
X = M. Denoting by G^(a|/?, t) the propagator of such a particle starting at a and ending, at time t, 
at /3, one has 


G+(a|/3,0 = E M 


e 2M^ 


Its LT wrt t, G~li{a\l3,s) reads 
G+(a|/3,s) = 

where we have used the identity 


k^l 


2sinh — inax(a,/?))] sinh /3)] 

sinh 

cos(fcx) 1 /-TTCOsh ((tt — x)a) 


+ k^ 2a 


sinh (tto) 


(46) 


(47) 


(48) 


As done before in (271, we introduce the two times tmax and r such that x(tniax) = M — e and 
x(t) = M — r. These two times break the interval into three sub-intervals (see Fig. |^. As shown in 
Fig. i two cases may arise: 0 <t < tmax or t > r > fmax- In this case, for the Brownian meander 
(BMe), these two configurations do not give the same contributions to the average DOS {pMe{r,t))- 
Using the same type of regularization procedure as used before (see Eq. ([2^), one has 


{pMe{r,t)) = lim 


1 


l-M-s 


s^O ZMe{s) 


dM / dU 


dxF 


(49) 


dr G^(e|M - r, T)Glf{M - r\M - e, tmax - T)Glf{M - £\xf: t - t^ax) 


-f f drG^(e|M-e,tmax)G^(M-£|M-r, r-ti„ax)G^(M-r|xF,t-r)j , 











14 



r —>■ 0 

r —>■ oc 

Brownian (BM) 

4r 

y TT r-^ 

Bridge/Excursion (BB/BE) 

4r 

4re“^’' 

Meander (BMe) 

4r 


Reflected Brownian 

4r 

8 /2 1 p-V 

Reflected Bridge 

4r 



Table 1 Asymptotic behaviors of average DOS for various constrained Brownian motions, both for small and 
large argument. 


where we have used the Markov property of BM and where the normalization constant ZMeis) is given 
by 


poo pt pM—e 

^Me(e) = / dM / dtinax / p GXj{Q\M - E G]^{M - e\x p ,t - 

Jo Jo Jo 


2e A / ^ j £ 0 ■ 


(50) 


To compute the numerator in (49), we take advantage of its convolution structure and perform its LT 


wrt t, using the expression of the LT of the propagator (47). After some manipulations, we obtain 
finally 


(pMe(r, t = l)) = 2v^ ~ ’ 

where = erfc(a;)/2 belongs to the family of functions studied in Appendix [ a} 


(51) 


3.4 Average DOS of constrained Brownian motions: summary and comparison 

It is useful to summarize and compare the results for the DOS for the BM and its variants which we have 
studied here, using a method based on propagators. In Fig. we have plotted (in lines) the results for 
the DOS for the Brownian motion ( |32[ ), Brownian bridge and Brownian excursion (42), the Brownian 
meander (51) as well as for the reflected Brownian motion and reflected Brownian bridge ( |154[ ), 
the study of which has been left in Appendix]^ We have computed numerically the DOS for these 
different constrained BMs - which have been simulated using the constructions described in detail in 
Appendix | e[ It is useful to remind that {p{r,t))/t = ,t = 1)) has the interpretation of a 

PDF, as dr(p(r, t))/t is the probability that the walker lies in the interval [xmax — r — dr, Xmax ~ x]- 
The main characteristics of these DOS are summarized in Tables[^and|^ A first interesting feature 
is that the small r behavior of the DOS is, at leading order, the same for all these constrained BMs, 
p{r, t = 1) ~ 4r. This indicates that the local vicinity of the maximum of the BM is insensitive to the 
boundary conditions in space (the presence of a wall at x = 0 either absorbing or reflecting does not 



numerics 


(r) 

^typ 

(r) 

rtyp 

Brownian (BM) 

= 0.7979... 

y TT 

0.5145... 

0.79 ±0.01 

0.52 ±0.02 

Bridge/Excursion (BB/BE) 

iyf = 0.6267... 

0.5 

0.62 ±0.01 

0.48 ±0.02 

Meander (BMe) 

i^/?(81og2 - 3) = 0.7975... 

0.4907... 

0.79 ±0.01 

0.48 ±0.02 

Reflected Brownian 

— ^0.7214... 

0.5212... 

0.72 ±0.01 

0.52 ±0.02 

Reflected Bridge 

iyf(41og2-l) = 0.5554... 

0.4907... 

0.55 ±0.01 

0.48 ±0.02 


Table 2 Ma in c haracteristics of the average DOS for various constrained Brownian motions: the average value 
(r) [see Eq. (@1 and the typical value rtyp which is the location of the peak of {p{r,t = 1)). The numerical 
result are obtained by averaging over 4.10^ realizations of independent RWs of 10^ steps. We refer the reader 
to Appendix E for a description of the algorithms which were used here. 
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affect it) as well as in time (the free BM and the BB, corresponding to periodic boundary conditions 
in the time direction, share the same local properties). On the other hand, for large argument, the 
DOS exhibits, in all the cases, a leading order Gaussian decay (see Table but with different rates. 
In particular, for bridges it decays much faster, oc , than for free processes, for which it decays as 

This also implies that the mean value 

(r) = / dr r {p{r,t = 1)) (52) 

Jo 

is larger for the bridge process than for the one which is free at the extremity of the time interval. 
Interestingly enough, although the mean value might differ notably from one process to another (see 
Table [^, the typical value r^yp, which is the location of the peak of {p{r,t = 1)), does not vary too 
much, with a typical value rtyp « 1/2 for all the processes (see Table and Fig. |^. 



r 


Fig. 7 (Color online) Plot of the average DOS for Brownian motion and its variants on the unit time interval. 
The symbols indicate the results of our numerical simulations (obtained by averaging over 10^ realizations of 


independent RWs of steDsl. The lines correspond to our exact analvtical results given in Eos. gl^lTl 
151 |l54 1. We refer the reader to Appendix [ e| for a description of the algorithms which were used here. 


4 Path integral approach 

As mentioned in the introduction, the average DOS is useful to compute the average of any functional 
of the maximum Oynax{t) = fgV(xmax — x(T))dT, see Eq. |2L However, if one is interested in 
higher moments of Ojnax(t), the knowledge of the average DOS is not enough. A useful approach 
to study the statistics of Omax(t), beyond the first moment, is to compute the Laplace transform of 
the full PDF of Oinax(t), namely {e~^Jo T(a:mi,x-a:(T))dT^^ from which the moments of Omaxit) can be 
obtained by successive derivation wrt A. In this section, we present a general approach, based on path 
integral techniques to compute this Laplace transform, for any function V(x). We develop this general 
framework both for the free BM and then for the Brownian bridge ~ which turns out to be easier to 
study, thanks to the Vervaat’s construction (see Fig. |^. We then apply these general methods to study 
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the full statistics of Pdos(7'j^) then to the study of the functionals Ta{t), corresponding to the case 
V{x) = x°‘, for the free BM and for the BB. Then, we focus on the special case a = —1, corresponding 
to Odlyzko’s algorithm, which, as we show below, leads to a quantum mechanical problem which is 
exactly solvable. 


4.1 General framework 
4 .. 1.1 Free Brownian motion 

To study analytically (exp[—A drl/(a;max ~ 3^('r))]), with A > 0, for an arbitrary function V{x), we 
first decompose the time interval [0,t] into two subintervals [0,fniax] and [tmax,i] where tmax is the 
time at which the maximum is reached. These two intervals [0,tniax] and [tmax,^] are statistically 
independent (as BM is Markovian), and the PDF of tmax is given by the arcsine law [JT] , 


P(imax) = 


fmax) 


(53) 


The process ?/(t) = a^max — x{t) is obviously a BM which stays positive on [0, t]. By reversing the time 
arrow in the interval [0, tmax] and taking fmax as the new origin of time, we see that 2 /(t) is built from 
two independent Brownian meanders (BMe): one of duration tmax and the other (independent) one of 
duration t — tmax (see Fig. [^. Therefore one has 


^ rdt,,ax<p(W)^(i-tmax) 
^0 

ip{T) = , 


(54) 

(55) 


where (•••)+ denotes an average over the trajectories of a BMe yir). In (551 the prefactor Xj^-KT 
comes from the PDF of t^ax (531. 


BM Two independent meanders 







. 


0 t ^max ^max t 


Fig. 8 (Color online) Illustration of the path decomposition used in the path integral method (54l: the BM 
is built from two independent realizations of a Brownian meander, one of duration tmax and the other one of 
duration t — tmax. 
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On the other hand the convolution structure in (54) suggests to compute its LT wrt t: 





(p{s) = / e 


(56) 

(57) 


The next step is to compute <^(s) in (57) using a path-integral method. For a Brownian meander, which 
is continuous both in space and time, this path integral method must be handled with care: as noticed 
repeatedly in the previous section, one can not impose simultaneously y(0) = 0 and y{0'^) > 0. This 
can be circumvented m by introducing a cut-off e > 0 such that ?/(0) = £ so that the Feynman-Kac 
formula reads: 


with 




^^^\e)dyF 

™°r(2/F|e-«ol£)dy^ ’ 


Hx 


2 dy^ 


+ XV{y) + Kvau(j/) 


(58) 


(59) 


where Kvau(y) is a hard-wall potential, Vwaii(y) = 0 for j/ > 0 and Fwaii(j/) = +oo for y < 0, which 
guarantees that the walker stays positive, as it should for a meander. We can easily compute the 
eigenfunctions of iJo = - 5^ + Kfaii(j/): 


i^k{y) 



sm{ky), 


(60) 


and 


{yrle = — I dk sin{kyF)sin{ke)e 


A straightforward computation yields 


dyFiyrle 


2e 


£—^0 


(61) 


(62) 


We denote by the wave functions of Hx associated to the energy (we will assume that the 
spectrum of Hx is discrete but the computation can easily be extended to a continuum spectrum): 





£—^0 2£ 



dyp J2'^niyF)'I^u{£)e 

n=0 


'' n=0 

00 

dyFY.'^:{yF)K{Q)e-^'^^ 

n—0 

dyp i9a;GT-(0, yp) , 

where we have used >F„(0) = 0, because of the absorbing wall in 0, and where we denote 

00 

Gr{x,y) = J2K{y)'^n{x)e-^^^ . 

n—0 



(63) 


( 64 ) 
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We can use these formulae (63 641 to calculate the Laplace transform ip{s) in (57) 

O pOO -I pOO 


(p{s) = / e ^'^ip(r)dt = 


/o 




dr e 


dyrd^^GriO^yp) = 


'0 


/o 


V2. 


0 


dyp Y. 


n—0 


= ^ J dyp d^Gs{0, yp) , 


KiyF)Kid) 

s + En 


(65) 


where Gsix,y) = 1+g" Laplace transform of GT-{x,y) wrt r. One thus recognizes 

that Gs{x,y) is the Green’s function satisfying 


[Hx + s] Gs(a;, y) = S{x - y) , 


( 66 ) 


such as Gs(a;, 0) = Gs{0,y) = 0. To compute Gs{x,y), we look for two complementary functions Us{y) 
and Vs{y) solution of the homogeneous equation [H\ + s] ipiy) = 0 with Ms( 0) = 0 and Vs{y —>■ oo) = 0. 
From Us{y) and Vs{y), we can compute Gs{x,y) as (42] 

2 

— Us{x)vsiy) iix<y , 

(67) 

2 

— Us{y)vs{x) iix>y, 

where 


Gs{x,y) = < 


W = u'^{x)vs{x) - Usix)v'^{x) 


( 68 ) 


is the Wronskian associated to u,(x) and v^(x), which is here independent of x. Finally we obtain 
from ( 1 ^ and ( 1 ^ 


POO JTs POO 

(p{s) = J e""^v5(T)dr = -^ J dyp u',{0)vs{yF) . 


(69) 


This formula (69) together with (56) and (57) allows one to compute (e a;(T))dT^ 


function V{x). To obtain explicit results from these general formulas, we need to analyze in more 


detail the Schrodinger operator in Eq. (59). This will be done, for some special cases, in section 4.2 


4-. 1.2 Brownian bridge 

In the case of a BB, the method presented above can be straightforwardly adapted to compute the 
Laplace transform of the PDF of functional of the maximum of the Brownian bridge, Omax.s {i)-, namely 
{e~^^o ^(^max,f 3 -ai 3 ('r))dT^^ principle, one could use the same reasoning as before, i.e. break the time 
interval [ 0 ,t] into [ 0 ,tinax] and [tmax,t] where tmax is the time at which the maximum is reached (see 
Fig. [^. The main difference is that, for the BB, the PDF of tmax is uniform P(tmax) = l/tj and 
not given by the arcsine law ( |53[ ) - this is a consequence of periodic boundary conditions in the time 
direction. There is however a simpler way to proceed, which makes use of the Vervaat’s construction 
(see Fig. [^, which allows us to map any functional of the maximum of a BB onto a (standard) 
functional of the Brownian excursion xe(t), from which path-integral techniques have already been 
developed liiTniiioi. Hence, generalizing the relation in (|35[) to more general functionals, we have 


f V{Xms^y,^B - XB{T))dT^= ( V{xe{t))At . 
7o Jo 


(70) 


Hence, using that identity in law (70), the LT of the PDF of any functional of the maximum of the 
BB, can be written as a path-integral [Rl lTT?ll2()| . Because we are dealing with a Brownian excursion, 
which prevents us to impose simultaneously xe(0) = 0 and a:£;( 0 ''') > 0 , this path integral method 
needs to be suitably regularized. This can be done, as explained before [see Eq. (39)], by using a 
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cutoff e such that a;£;( 0 ) = xsit) 
a limiting procedure: 


e and computing the statistics of these observables (701 by using 




e^o {e\e~^»*\£) ’ 


(71) 


where H\ is the Schrodinger operator defined above (59l. Note that the denominator in (711 does not 
depend on the functional at hand, i.e. it is independent on V(x), and it is readily computed from (61 \ 
to be 


(e|e ~ , e ^ 0 . (72) 

This formula 0 is quite general and it allows us to compute the statistics of a wide class of functionals 
of the maximum of the Brownian bridge, as we illustrate it below. 


4.2 Applications to some specific functionals of the maximum 

In this section, we present some concrete applications of the above path integral formalism, both for the 
free BM as well as for the BB. We first illustrate the method on the computation of the full statistics of 
the DOS, providing a detailed derivation of the results recently announced in Ref. [^. Then, we study 
the special family of functionals corresponding to V(x) = a;“, whose applications where discussed in 
the introduction. Finally, we provide a full detailed study of the special case a = —1, i.e. V{x) = 1/x, 
which corresponds to the analysis of Odlyzko’s algorithm. 


4-.2.1 Full statistics of the DOS 


The main body of results for the statistics of the DOS [see e.g., Eqs. ( |8^ and (83) below] were recently 
announced by us in a short Letter [23]. In this section we provide a detailed derivation of these results 
using the path-integral framework presented above. 

The case of free BM. To compute the statistics of p(r, t) = S(xmax — x{t) — r)dT, we apply 
the above formalism (56 57 69) to the special case V{x) = S{x — r). In this case, the two independent 
solutions are Us{x),vff^ of [Hx + s]'!/’(a;) = 0, with Hx = —(l/2)d^/da:^ -I- XS{x — r) + t4ran(a^) which 
reads simply here, for x r: 


— ^if"{x) + s tfix) = 0 , for X G [0, r[U]r, -|-oo[ 


with the following boundary conditions: 

zis(O) = 0 

Usir+) = Us{r~) 

i«(r+)-<(r-)) =Xu,{r) 


lim Vs{y) = 0 

y—yoo 

Vs{r+) = Vs{r~) 

i«(r+) -u((r-)) = Xv,{r) , 


(73) 


(74) 


where the boundary conditions in r result from the presence of the delta peak in Hx = —\df/dx^ + 
X5{x — r). Hence, Us(x) and Vs(x) are given by 


f Asinh(v^x) 

Us(x) I A ^sinh(-\/2sx) -I- sinh(-\/2s(a; — r)) j 

() - f ^ |'g-V2l(x-r-) ^ (x/^(r - x))) X < r 


X <r 
X > r , 


(75) 

(76) 


X > r, 

where A and B are normalization constants which are irrelevant here. The Wronskian is thus given by 


W = u'g{x)vs{x) — Us{x)v'g{x) = ~ ® j A-f 


(77) 
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With these formulas ( [^ [M| fnl , we finally find 

+ X (l - 

y2s + A (^1 - 6-2^25’-j 


(p{s) = 


(78) 


which, combined with Eq. (56), yields finally the formula : 




1 f v^+A(l-e-^0^ ' 

V2s + A(1 - 


(79) 


From this expression (79), one can obtain the moments of arbitrary order (p^(r, t)), for fc € N by 


expanding it in powers of A. It yields: 

pOO oo 

/ dte-"‘(e-^^(’’’*)) = l + ^(-A) 


fc = l 


Q -v/^r 

fe»^^_(l_g-V27r-)2(i_^^g-V27r)(i_g-2V27r)/c-2 ^ (gQ) 
(2s) 2+1 


This expansion (80) in powers of A yields the LT of the moments {p^{r,t)) wrt t as: 

r(p'=(r,t))e-«‘dt = ^(_i)'(fcyi) ((2Z - i)e-(2'+i)v^’' + (k - 2(1 + l))e-^('+i)-^’').(81) 

Jo (2s) 2 +^ V / 


(2s) 

It is then possible to invert this LT using the functions presented in Appendix [a| to obtain 

fc-i 

(p'=(r,t = 1)) = 8 k\ Y.{-iy{’^-%2l + !)<?('=+')((2Z + l)r) + (fc - 2 {l + 1))^^’^+^\2{1 + l)r)] .(82) 


1=0 


For k = 1, this yields back the result obtained in (32). 

By inverting the LT wrt A in (79), we obtain the full PDF Pt{p, r) of the DOS (Q, as a function of 
p, for different values of the parameter r. 


e '*‘Pt(/ 0 ,r)dt = (5(p) 


(e 


— y/^r 


- 1 )^ 


pv^eV^y 

g 2=inh(v^p) 


e 2 


pe 




s(l + e-'/^'-)2 cos^a 4sinh (rV^) sinh 


.(83) 


It has an unusual form with a peak oc 6 {p) at p — 0, in addition to a non trivial continuous background 
density pt{p,r) for p> 0. Hence one has 

Ptip,r) = Fw{r,t)5{p) +pt{p,r ), (84) 

where Fy/{r,t) = Prob.[IT(t) < r], given below in ( |^ , is the probability that the width W{t) = 
max^g[o^t] x{t) — min,-g[o_t] x{t) is smaller than r. This can be understood because if W{t) is smaller 
than r, the amount of time spent by the process at a distance within [r, r + dr] from the maximum is 
0 (see Fig. J^, yielding the delta peak at p = 0. Indeed we can check that the coefficient of the term 
cx 5{p) in (83) is the LT wrt t of 


Fw{r, t) = 1 + 4^(—l)*erfc(P/-\/^), 


(85) 


1=1 


which corresponds precisel y to the distribution of the width (or the span) of BM |43[l44] . 

On the other hand, in (84), pt{p,r) = p\{pj\fi,rj^fi)j\ft is a regular function of p, for r > 0 (see 
Fig. [|), and has a more complicated structure. We obtain an explicit expression of its LT wrt t given 
by the second term of Eq. (83). In Fig. |^we show the results of pi{p,r) obtained from numerical 
simulations (averages are peirormed over 10^' samples) for three different values of r. We see that 
they are in perfect agreement with our exact formula (83). The analysis of the distribution pi(p, r). 
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Fig. 9 (Color online) Plot of Pi(p, r) as a function of p for different values of r. The solid lines for p = 0 
represent the 5{p) peak, oc 5(p) in Eq. ( |83[ ). The dotted lines correspond to our exact analytical results for the 
background density pi (p, r) in (831 - where the inverse LT wrt s has been performed numerically (in green for 
r = 2.5, purple for r = 1.2 and blue for r = 0.2) - while the red dots Indicate the results of simulations. On 
the 2 = 0 plane, we have plotted the exact mean DOS in Eq. (32l. 


including its asymptotic behaviors when p —>■ 0 and p —^ cx) was carried out in Ref. |29| and we refer 
the interested reader to the supplementary material of Ref. |29] for more details. 

The case of a Brownian bridge. In this case, we remind that, as it can be seen using the 
Vervaat’s construction, the DOS for the BB coincides with the DOS for the Brownian excursion, which 


in turn, coincides also with the local time of the excursion [see the identity in law in Eq. (35l]. The 
full PDF of the can be obtained from the general formula (0, applied to the special case 

V{x) = 6{x — r): 


(e-Vn(r.i)) ^ 


e e 


.-H,t 


e) 


1 d2 




( 86 ) 


where we recall that Vwan(2() is a hard-wall potential, Twan(2() = 0 for a; > 0 and = -|-(X) for 

a: < 0 - hence Hx is defined on [0,-|-oo[. The numerator in (861 can be computed by expanding the 
matrix element on the eigenfunctions tpk {x) of Hx which are given by 


0 a: < 0 

Hxipk{x) = y sin(A;a;) 0 < a: < r 

a{k, A) (sin(A:a:) -b ^ sin(fcr) sin(fc(x — r))) r < x, 


(87) 


where the amplitude a(fc. A) ensures the normalization condition: 


\a{k,X)\-^ = - 


1 -b — sin(A:r) ( cos(fcr) -b — sin(fcr) 
k \ k 


( 88 ) 


We finally obtain 





2t3 

TT 


1 -b ^ sin (kr) sin (fcr) -b cos (fcr)) 


'dk 


(89) 


Note that from the identity in law in (35) this expression (|89|) yields also an interesting relation for 


the local time of the Brownian excursion, which we have not seen in the literature. Furthermore, by 
expanding this formula (89) in powers of A, it is possible to obtain, using yet another method, the 
moments ([pb(c, f)]^) of arbitrary order k gN and recover the expression given above (42). 
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On the other hand, by studying the large A behavior of (891, which is of order 0(A°), we can show 


that the PDF of pB(r, t) has an expression similar to, albeit different from, the one for BM in ( |M| : 

Pt,B {p, r) = Fw,b {r, t)S{p) +pt,B{p,r) , (90) 

where distribution function of the width of the BB |45| : 


Fw,b (r, t) = 1 + 2 ^ < 


1=1 


1 - 


t 


(91) 


while Pt,B{p,i^) is now a different distribution. In particular, one can show that (421 yields back the 
complicated though explicit formula for pt^B{p,f) found in Ref. |39LI46| using a completely different 
method: 


OO j-l / . V 

Hk+2{p + ‘^i'j)/kl 

1 = 1 fc=o '' 


j=l k=0 

where the iF„’s are Hermite polynomials defined by 

[n/2\ 


Hnix) = n\ Y, ^ 


(- 1 )* 


i=0 


2U\{n-2i)\ ’ 


where [xj denotes the largest integer not larger than x. We refer the interested reader to Ref. 
more details on the distribution pt^Bip,x), including its asymptotic behaviors. 

4-2.2 The case V{x) = x°^: first and second moments of the functionals 


(92) 

(93) 
3 for 


We now apply this path integral formalism to the functionals Taft), for free BM (12l, and Ta^Bft) the 


associated observable for the BB ( |13| ). We showed previously that the first moment can be obtained 
directly from the corresponding average DOS. This simple method can not be easily adapted to compute 
higher moments of Ta{t) or Ta^Bft). We show here how to compute these moments using the path 
integral formalism developed above. We will treat separately the case of free BM and BB. 

The case of free BM. It is convenient to start from the following formula, obtained from the 
combination of Eqs. ([56][^|6§: 


nOQ 1 / nOO \ 2 

= _ M dy5,G3(0,j/)j , 

where we recall that Gs{x,y) is the Green’s function of the following Schrodinger equation 


~dx^ + AR(x) + Vwaii( 2 :) + s 


Gs{x,y) = S{x - y). 


(94) 


(95) 


Of course, there exist very few instances of potential V(x) for which this Eq. (95) can be solved 
exactly (see below). However, the first moments of Taft) can be extracted, in principle, for a generic 
V{x). Indeed, from (94), we see that the moments of Ta{t) are obtained by successive derivations 
of the right hand side of (94) wrt A, evaluated in A = 0. These successive derivatives can in turn 
be expressed as combinations of the successive derivatives of Gs{x,y) wrt A. Indeed, if we write the 
following expansion: 


A” 


Gsix,y) = V —rfn{x,y) , 

n\ 


(96) 


n—0 


one has from (94): 


kl 


{Tt{t))e-^*dt =i-l)Y Y^^r^bk- 


n 7 I 


n=0 


1 

n! 


dxfn{x,y)\^^^ dy 


(97) 
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Hence, we need to solve perturbatively the equation for the Green’s function (951 in powers of A to 
compute the functions fnix,y) in (96l. First it is easy to obtain fQ(x,y) as 


foix, y) = _ e-V 2 A(^+ 2 /)) _ 

V 25 


(98) 


Furthermore, one can show that the functions fn{x,y) satisfy the following recursion relation: 

nOO 

fn{x,y) = -n dzfo{x,z)V{z)fr,-i{z,y) , n>l. (99) 

Jo 

The recursion relation ( |M| can be solved formally in the closed form 

n .oo "-1 

fn{x,y) = {-irnlYl / <iZifo{x,Zi)Y[V{zi)fo{Zi,Zi+i)V{Zn)fo{Zn,y)- (100) 


i=l 


Hence, from Eq. (971 together with (1001 one can compute the moments of arbitrary order of any 
functional of the maximum of the free BM. Note that this technique is a generalization of the method 
based on propagators which we used before for the computation of the average DOS, corresponding to 
a particular potential V{x) = S(x — r). Here we apply this formalism to the cas e Vj x) = x°‘. 

We first compute the first moment, using Eq. (971 for k = 1, together with (lOOl. We obtain: 


pOO / pOO \ / pOO 

J dte"®*(T„(t)) = - dyd^fi{x,y)\^^^j dyd^fo{x,y)\ 


x—0 


( 101 ) 


Using the expressions of fo{x,y) in ( [M| and of fi{x,y), obtained from Eq. ( 100[ ) together with the 
fact that dxfo{0,y) = we obtain (performing the change of variable z —)■ \^z): 


dydxfi{0,y) = - dy dzdxfoi0,z)z°‘foiz,y) 
Jo Jo 
2 


2 


Q!+3 


pOO pOO 

dy dze-"z“(e-l"-^l-e-("+^)) 
Jo Jo 


( 2 s)- 


-r(a + l)(2-2-“). 


Hence, combining (101 1021 and using that dydxfo{x,y)\ = a/2/s we arrive at 


4 

/ dte-«*(T„(t)) = —^T(a + l)( 2 - 2 -“), 

Jo {2s)^ 


which after Laplace inversion, yields 


(Tcit)) = 


( 2 t)i+t( 2 - 2 -“)T(l^) 

(2 + a)^ 


( 102 ) 


(103) 


(104) 


This formula coincides, as it should, with the formula obtained directly from the average DOS in (33l. 
We now compute the second moment, using Eq. (97) for k = 2: 

poo / poo \ 2 / poo \ / poo \ 

Jo ^ dydxfii.0,y)j + <^ydxfoi0,y)j\^J^ dydxf2i.0,y)j . (105) 

We then have, using (98 1001, and performing the change of variables with the substitution u —> zi , 

V —>■ Z2- 

pOO pOO pOO pOO 

/ dydxf2{0,y) = 2 dy dzi dz2 d^foiO, zi)z^ fo{zi, Z2)z2 fo{z2,y) 

Jo Jo Jo Jo 

A pOO pOO pOO 

/ dy du due-“u“(e-l“-’’l-e-(“+’’))u“(e-l’’-^l-e-(''+^)) 

Jo Jo Jo 




5+2a 


•\/ 2 s 


5+2a 


poo poo 

Ca,withC'„=/ du due-"M“(e-l“-''l-e-(“+“))u“(l-e-’’) (106) 

Jo Jo 
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where Ca can be explicitly computed as 

r{2 + 2a + n) 




^ n!22+2“+"-(l + a + n) 23+2“(l + a) 




‘2^‘2-\-‘2o 


Finally, using (98 102 105 1061, we obtain 


4 

I dte-^\T^{t)) = {r{a + 1)^(2 - 2-“)2 + 4C„) , 


which after Laplace inversion yields: 

4.2 + ' 

imt)) = 


23“r(3 + a) 


F(a + 1 ) 2 ( 2 “ - l)( 2 “+i - 1 ) + 

F(2 + 2a + n) \ 
^n!2i+”(l + a + n) ) 


r(3 + 2a)(4“+i - 1) 
4(l + a)2 


+E 


(107) 


(108) 


(109) 


The last sum over n can be finally expressed in terms of incomplete beta function, which yields the 
formula ( |16[ ). In particular, by taking (carefully) the limit a —> — 1 in the above expression (109), we 
recover the result of |33] 


lim (r2(t = l)) = ^+41og(2)2, 
a^-l 3 

which was obtained by the authors of [33] using a completely different method. 


( 110 ) 


Brownian Motion (BM) 



Brownian Bridge (BB) 



a 


Fig. 10 Left panel: Plot of {T^{t = 1)), as a function of a, for the BM, as given in Eq. (1091. Right panel: 
Plot of {T^^g{t = 1)), as a function of a, for the BB, as given in Eq. (1211. In both panels, the dashed black 
lines indicate the asymptotic behaviors discussed in the introduction. In particular, for q = — 1 for the BM 

(corresponding to Odlyzko’s algorithm), we recover (T^__i(t = 1)) = ^ + 41og(2)2 [33] . For a = 1 for the 
BB (corresponding to the Airy distribution), we recover = 1)) = ]^ |19II20| . 
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The case of BB. In this case, the equivalent of Eq. (941 is the formula derived in Eqs. (711 and 


(721 which reads 




Taking the LT of the above relation wrt t, one obtains: 


M)) = ^d,dyG,{x,y) 


x—y—0 


( 111 ) 


( 112 ) 


where Gs{x,y) is the Green’s function of the operator in (95l. From (112), one deduces the following 
relation 


e-"*r3/2/ 


/ ^max,B - XB{r))dT 

Jo 




X—y—0 


(113) 


where fk{x,y) is given in Eq. (100). From (113), specified to the case V(x) = x°‘ and fc = 1, we obtain 
the first moment {Ta^sit)) as 


/*oo j o2 


x—y—0 


Hence, using (98 100) together with the change of variable z —>■ V^z, we obtain 

pOO j ^OO 

j J dzd^fo{ 0 ,z)z°‘dyfo{z, 0 ) 


(114) 


(115) 


TT 4 


v2s -'0 

which, after Laplace inversion, yields immediately 

^l+a/2 


dze-^^z“ = 


TT 2^ 




.T(l + a), (116) 


(117) 


which coincides with the formula obtained above from the average D OS ([44| . 

We can also compute the second moment (Tj ^(t)} by using Eq. (112) for the case k = 2, yielding 


dte = ^d^dyf 2 { 0 , 0 ) 


(118) 


where f 2 {x,y) can be obtained from (100) with V{x) = x°‘. Using (98 100), we obtain, performing 
the changes of variables u —>■ v —>■ •\/^Z 2 , and after some manipulations 

/'OO pOO pOO 

/ dte~'‘H~^^‘^{T^g{t)) = V^ dzi dz2dxfo{.0,zi)z^fo{zi,Z2)z2dyfo{z2,0) (119) 

Jo ’ Jo Jo 


( 2 r {2 + 2 a) 


22«y^'“+H 1 + ' 

which, after inverse Laplace transform, yields 


-U(l + a)' 


/ 2r(2 + 2a) 


( 120 ) 


( 121 ) 


In particular, we can check that {T^i = 1)) = coincide with the known result for the area 
under a Brownian excursion, i.e.the second moment of the Airy distribution imiin]. 
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4-2.3 The exactly solvable case V{x) cx 1/x and application to the Odlyzko’s algorithm 


In this section, we apply our general formalism to the special case V{x) = 1/x, which is relevant to 
analyze the distribution of the cost of the optimal algorithm to find the maximum of a random walk 
(see Eq. (|4|), as shown in Refs. |301I33| . We apply this formalism separately both to the free Brownian 
motion and then to the Brownian bridge. 

The case of free BM. In this case, we want to construct the Green’s function Gs{x,y) in (671 
which can be built from two independent solutions Us{x) and Vs{x) of the following Schrodinger equa¬ 
tion 


1 d^ A 

2 dx^ X 


'ip{x) = 0 , s > 0 


( 122 ) 


with the asymptotic behaviors 


,(0) = 0 & Vs{y —>■ oo) = 0 . 


(123) 


It turns out that the above Schrodinger equation (122) can be solved in terms of hypergeometric 
functions. The solutions Us{x) and Vs{x) satisfying (1231 read 


Us{x) = A e iFi ^12, 2-\/^a;^ , 

Vs{x)=Be-'^^^xu(^l + ^,2,2V^sx'^ , 


(124) 

(125) 


where iFi(a,b,x) and U{a,b,x) are confluent hypergeometric functions (respectively Rummer’s and 
Tricomi’s hypergeometric function) and where A and B are two irrelevant numerical constants (note 
that Us{x) and Vs{x) are also known as Whittaker functions). The Wronskian W reads 


W = u'^{x)vsix) - Us{x)v'^{x) = 


AB 


2AT(A/v^) ■ 


(126) 


The function in Eq. (691 reads in this special case (124 1251: 


Note that one can check from (1241 that Us(0) = 0, as i t should [see Eq. (123)] and that u(,(0) = A. 

(127) 


^{s) = 2^/2 AT 


(«)f 


e-^^y yU { 1 + ^^,2,2V^sy] dy . 






After some manipulations, the integral over y in Eq. (127) can be evaluated as (see Appendix D.l) 


1 


VS \'j2s J 


A 


OO OO / V 

G{x) = 2 ^(-l)'=+iC(fc)cr'= , with C(fc) = (1 - 2-i-'=)C(A) = E M 


n+1 

k 


(128) 


fe =0 


Hence, from the general formulae (56 57), together with (128), one obtains the explicit formula for 
the case V(a:) = 1/x: 


3-«t(e-A/o‘ 


4 A" ~ - 

E ■ 

n=0 ’ /c=0 


(129) 


Expanding the left hand side of Eq. (129) in powers of A we obtain the moments of Ta=-i{t = 1) of 
arbitrary order, as announced in the introduction (18), 


{T^=.At = i)) = r 


fc + l\ 2t+2 


E C{m)C{k-m) . 


m—0 


(130) 
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This allows us to recover in a completely different manner the results obtained in |33| by probabilistic 
tools. For completeness, we mention that the authors of |33| obtained an explicit expression of the 
PDF p{s) of Ta=-i{t = 1) as (see Theorem 4.2 of that paper) 


^ = J - y)dy , 

where the functions 9i{x) and 02 {x) are given by 

"■W'l:/ 


a 


92ix) = ^ X] “ n^n‘^x)e~ 


(131) 


(132) 


(133) 


The case of the BB. In this case, the starting point of our analysis is the general formula given 
in 


exp 


-A 


dr 


= lim 


{e\e ^^*|e) 


- xb{t))\ / E^o (e|e--f^o*|e) 
where Hx is the Schrodinger operator defined in (|59[), which reads here 


~ ~2~^ ^ K.aii(a;) , A > 0 , 


(134) 


(135) 


whil e the small e behavior of the denominator of (134) is given in (72). To compute the numerator in 
Eq. (134), we expand the matrix element on the eigenfunctions \4'e) if Hx, which satisfy 


Hx\(I>e) = E\4>e) 


(136) 


where the eige nvalu es E > 0 form a continuous spectrum, as there are no bound states here [we recall 
that A > 0 in (135)] and where the eigenvectors satisfy the boundary condition 


lim 4>e{x) = lim {x\4)e) = 0 . 

x—>0+ a:—>0+ 


(137) 


The general solution of the Schrodinger equation (135 136) reads 

(t>E{x) = CExe-^^^ iFi - i-^,2,2ixy/^^ + (Iex - i-^,2,2ixV^^ ,(138) 

where we recall that iFi{a,b,x) and U{a,b,x) are confluent hypergeometric functions. The boundary 
condition at a: = 0 (137) imposes that d-E = 0. The remaining task is to compute the normalization 
constant ce such that 


noo 

/ (f>*E,{x)(f>Eix)dx = S{E — E') . 
Jo 


Using a formula given in Landau-Lifshitz (see Appendix D.2), one can show that 

4A 


|c£;P = 




(139) 


(140) 




Therefore, Eq. (134), together with Eq. (140) yields 

pt 


exp 


-A 


dr 


0 (^max,B 


= 4,/^e/2^ 


^-Et 


0 (exp [27tX/{V^)] - l) 


d£;. 


(141) 


It is then possible, using residue theorem, to invert the LT wrt A in ( |141| ) and then perform the 
remaining integral over E. This finally yields the PDF of Ta=-i = 1) given in Eq. ( p0| ). The 
expression for the moments given in Eq. (19) then follows straightforwardly. 
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5 Conclusion 

In this paper, we have presented several analytical tools to study functionals of the Brownian motion 
and its variants, in particular Brownian bridge and Brownian excursion. These tools include (i) a 
“paths counting” method, relying on propagators of BM with appropriate boundary conditions and 
(ii) a suitably adapted path-integral method, which allows us to recast the study of functionals of the 
maximum of BM into the study of quantum mechanical problems. The first method (i) is conceptually 
quite simple and allows us to obtain in a rather simple manner the mean value of any functional of 
the BM, while the second method (ii) is better adapted to compute the full PDF of such functionals. 
We have used these methods to calculate the statistics of the density of near-extremes, or density of 
states (DOS), for Brownian motion p{r,t) and its variants. In particular, from the mean DOS {p{r,t)), 
one can compute the average value of any functional of the maximum of BM. Then, we provided 
a thorough study of functionals of the form Ta{t) = fg(xmax — a;(r))“dr, with a s] — 2,-|-oo[. As 
a is varied, Ta(t) interpolates between various physical observables, as discussed above. We have 
obtained an exact expression for the two first moments {Ta{t)) and {T^{t)) both of which exhibit a 
non-trivial, non-monotonic behavior as a function of a. Thanks to the path-integral method, when 
the associated quantum problem can be solved exactly, it is possible to obtain an explicit expression 
of the Laplace transform for the full PDF of Ta{t) or from which moments of arbitrary 

order and in some cases the full PDF can be computed. We have worked out in detail the case 
a = — 1, corresponding to V{x) (X \/x, which corresponds to the cost of the optimal algorithm (due 
to Odlyzko’s) to find the maximum of a discrete RW. In this case, we provided an explicit expression 
for the moments of arbitrary order , recovering by physical methods the results obtained 

in Ref. |88| by completely different probabilistic approaches. Furthermore, we have generalized these 
results to functionals of the Brownian bridge, Ta^sit). In particular, we argued that, for a = —1, the 
random variable describes the cost of the optimal algorithm (i.e. Odlyzko’s algorithm) for 

the search of the maximum of a RW in a bridge configuration and computed explicitly its PDF as well 
as its moments of arbitrary order. 

Several interesting questions are left open. For instance, here we have studied the case of a single 
Brownian motion and it would be interesting to extend this study to the case of multi-particle systems, 
where there are > 1 walkers, which could be independent or instead interacting, as in the case of non¬ 
intersecting (vicious) walkers, whose extreme value statistics have recently attracted some attention 
m- Finally, we have treated the case of Brownian motion, which is continuous both in space and 
time, and it would be interesting to extend these results to random walks, which are discrete in time. 
This would, in particular, allow us to study the DOS of Levy flights, whose behavior is expected to be 
qualitatively different from Brownian motion. 


Acknowledgements We acknowledge support by the Indo-French Centre for the Promotion of Advanced 
Research under Project 4604 — 3. We acknowledge a useful correspondence with Philippe Chassaing. 
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A Some useful functions 


We introduce the family of functions j G N, which satisfy 




f 

JQ 




*dt. 


These functions can be obtained explicitly by induction, using |48| 


(142) 


V 27r 


roo 

{x) = / $^^\u)Au . 

J X 


The first functions can easily be computed as 





<?«(x) = ierfc(^^) , 

^ _ — -xerfc • 

More generally, one can show |48| that they can be written in the form 

+'7t(a;)erfc(^) , 


(143) 

(144) 

(145) 

(146) 

(147) 


where Pj{x) and qj{x) are rat iona l polynomials of degree j — 2 and j — 1, respectively, for j > 2 |48| . We refer 
the interested reader to Ref. [48| for efficient algorithms, which can be implemented numerically, to compute 
these polynomials in a systematic way. 


B Average DOS for reflected Brownian motion 


Using the method based on propagators, explained in section 3.1 see Eq. (271, we can also compute the 


average DOS for the reflected Brow nian motion xr{t) which is tire absolute value of the Brownian motion. 


xr{t) = | 2 :(r)|. The expression in i27[, see also Fig. 
of the reflected Brownian motion such that xr{t) < 


indicates that we need to compute the propagator 
or equivalently —M < x{t) < M. Therefore, we 


compute the propagator of a Brownian particle confined in a given interval [—M, M] with absorbing boundary 
conditions both in x = —M and x = M. Denoting by G^j{a\l5^t) the propagator of such a particle starting at 
a and ending, at time t, at jS, its LT wrt t is given by 


2 sinh I 

^V^(M — max(a, /3)) 

) sinh I 

+ min(a, /3))^ 

sinh ( 

'•JTs2M^ 

1 


(148) 


In order to compute the average DOS {pR{r,t)) for the reflected BM, we evaluate the “number” of Brownian 
trajectories satisfying the following constraints: the process reaches its maximum M or its minimum —M at 
time text, passes through M — r or —M + r at time t and end in xr € [—M, M] at time t. The total number 
of such trajectories is then obtained by integrating over xf,M and text. When dividing the time interval [0, t] 
into three parts delimited by r and text, 8 different cases may arise: r < text or r > text, x(text) = ±M and 
x{t) = ±(M — r). Using the invariance of the process under the reflection symmetry x —1 —x we have to 
consider only four different cases (each one with a multiplicity of 2): 


2 poo nt pM 

{pR{r,t)) = lim / dM / dtext / dxF 

e^0Zr{€) Jq J-M 

[ J “ drG^(0|M-r,T)G^(M-r|M-e,text-r)G£(M-£|xF,t-text) 

+ r drGf_f(0|M-£,text)G£(M-e|M-r,r-text)Gff(M-r|xF,t-r) 

''text 

ptext 

+ / drGfKOk - M, T)G^(r-M|M-£, text-r)Gf_f(A7-£|®F,t-text) 

Jo 

+ r drGf,(0|M-£,text)G£(M-£|r-M,r-text)Gff(r-M|xF,t-r)l 

Jt^^^ -I 


(149) 
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where we have used the Markov property of BM and where Znit) is the normalization constant (such that 
/o°°dr(pH(r,t)) = t) 


POO pt pM 

Zr{£)=2 / AM / dtext / AxfGm(A\M - £,U^t)GM{M - e\XF,t - tey,t) ■ 
Jo Jo J-M 


(150) 


The normaliz atio n is easily computed as Zii{e) ~ 2£^, when e —>■ 0. Using th e sam e kind of calculations 
as in section 3.1 - exploiting the convolution structure of the integrals in Eq. (1491 - we find, after some 
manipulations 


{pR{r,t = 1)) =8^ 


(- 1 ) 


n + 1 


3 - 2n + 12n2 + 8n3 , 

n=0 \k=0 


^(3 + (-l)''(2n + fe)''+^)(2n + kf<P^'^\{2n + k)r) ) , (151) 


where 4>^^\x) is given in Eq. (146l. 

Similarly, we can study the UUS of the reflected Brownian bridge xrb{t), which is the absolute value of a 
Brownian bridge xr b(t) = |a;sfl(T)|. The calculation of the DOS in this case is very similar to the case of the 
free reflected BM in (1491 without the integral over xf which is set to xf = 0. Using time reversal symmetry, 
we can show that the average DOS {pRB{r,t)} is given by 

4 /■* /■* 

{pRB{r,t)) = lim— / AM / dtext / dr 
e^oZ Jq 

[ G£(0|M - e,text)G£(M - e|r - M,r - U^t)Gl,{r - M|0,t - r) 

+ G£(0|M - e,text)G£(M - e|M - r,r - text)G£(M - r|0,t - r)] 
where ZRB{e) is the normalization constant, given by 

POO pt 

I AM 
Jo Jo 


Zrb{s) = 2 


dtext Gm{0\M — e, text)G^(Tf — e|0, t — text) • 


(152) 


(153) 


The normalization is easily computed as ZRBie) ~ 2e^/(VSrrt), as e —>■ 0 and eventually the average DOS 
{pRB{r,t)) is obtained as: 


( OO 

4 n(—1)"^^^*-^^ (2nr) — (2r) 

n = 0 


(154) 


where $^^\x) is given in Eq. (145l. 


C Odlyzko’s algorithm 

C.l Main ideas behind Odlyzko’s algorithm 

To get familiar with this algorithm, it is useful to consider a simpler search algorithm, denoted by u, belonging 
to An (that denotes the ensemble of the algorithms that find the maximum of a random walk of n steps), 
which proceeds as follows: u probes always the random walk at the step where the upper envelope of the (still) 
possible trajectories reaches its maximum. This algorithm u is based on the idea that, as illustrated in Fig. 
|11| if Xm and X^+k have been probed, then the searcher knows for sure that, between step m and step m + k, 
the position of the random walker can not exceed {Xm + Xm+k + k)/2. This can be shown as follows. Let us 
denote by n+ the number of up-steps (-1-1) and n_ the number of down-steps (—1) between step m and step 
m + k. Then n+ and n_ satisfy the equations 



n+ + n_ = k 

71+ Tl— = ^m+k • 

(155) 

(156) 

Hence one has 

^ ^m + k Xm “h k 

(157) 


Xm Xm-\-k k 

(158) 

Therefore the position of the random walker can not exceed Xm + n+ = {Xm + Xm+k + k)/2, as shown in 
Fig. This simple algorithm is illustrated in Fig. [^on a realization of the RW for n — lA steps. This basic 
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Fig. 11 Illustration of the main idea of Odlyzko^s optimal algorithm. The R.W can not exceed T T 

k)/2 between m and m + k. If this quantity is smaller than a new probe between m and m + fc is useless. 



Fig. 12 An example of the algorithm u for finding M 14 = 7 for a RW in 4 probes, (a) Typical realization of 
a 14 steps RW, for which we want to find the maximum. We know without any probe that 0 < M 14 < 14, 
and if M 14 = 14, the maximum would be at position 14 (RW with only +1 jumps) so we probe the position 
14. (b) The first probe shows X 14 = 6 and we know now (see Fig. 11 1 that 6 < M 14 < 10, and if M 14 = 10, 
the maximum would be at position 10 (dashed line) so we probe the position 10. (c) The second probe shows 
Xio = 4 and we know now that 6 < M 14 < 7, and if M 14 = 7, the maximum would be at position 6 or 
13 (dashed line) so we probe the position 6 . (d) The third probe shows Xg = 4 and we know now that 
6 < Mi 4 < 7, and if M 14 = 7, the maximum would be at position 13 (dashed line) so we probe position 13 and 
find the maximum M 14 = X 13 = 7 in 4 probes. 


idea is at the heart of the algorithm proposed by Odlyzko. 

Here we also want to explain briefly the occurrence of this particular fu nction al of the maximum 1 in Q, 
in the analysis of this optimal algorithm, following the line of reasoning of mm- To understand this, let us 
consider a traveler, moving on a line, its position being denoted by y. Suppose that its velocity v{y) at position 
y is bounded by some function z{y), such that 0 < v{y) < z{y). Then the time t to reach the point x starting 
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from the origin satisfies the bound 


. r ^ > r ^ 

Jo ^(y) ~ Jo z{y) ■ 


(159) 


Now let us consider an algorithm a, its cost being C{a) and denote by mi,... ,mc{a) the steps at which the 
RW has been probed by the searcher - which has eventually found the maximum after C{a) probes. To 
be sure that the maximum is not in the interval [rrii, rrii+i], the potential maximum of the RW between these 
two steps, which is {Xmi + Xrm^i + rrii+i — mi)l2 (see Fig. Illll, must be smaller than M„ (by definition of 
the maximum). Hence this yields the following inequality 

mi+l TTli ^ 2Alfi XfYi^ . (^^9) 


Not ice that rrn+i — rrii can be seen as the velocity v{mi) of the alg orith m at point rrii. One can further argue 
|30| . using the fact most of the RWs are “slowly varying” [see Eq. (1631 below], that 2M„ — Xm- — Xm^^^ ~ 
2{M„ — Xrrii) when n is large. Hence 


Zk — 2(M„ — Xk) 


(161) 


can be viewed as the speed limit at step k of the random walk. Finally, by analogy with (1591, C{a) satisfies 


n - 1 


2 Mn — Xk 


(162) 


which in the continuum limit yields the functional of the maximum 7 in Eq. (j^. It is rather clear that these 
heuristic arguments leading to Eq. (162 I can be straightforwardly extended t^he case of the Random Walk 
bridge, Xi^B, which is a RW conditioned to start and end at the origin Xo,b = Xn,B = 0. Of course in this case 
the maximum M„ in (1621 is then replaced by the maximum of the Brownian bridge = maxi<i<jv Xi^B- 


C.2 Description of the Odlyzko’s algorithm 

Here we describe in more detail Odlyzko’s algorithm which finds the maximum of a random walk Xi+i = Wi 1 
with equal probability 1/2 (starting from Xq = 0). Let c be a positive real number, which is sufficiently large. 
The algorithm is essentially based on the fact that most of the RWs has “slow variations” (SV), i.e., check the 
identity [30| : 


\Xi+k — Xi\ < cy/k log n , Vi, k with i + k < n . 


(163) 


Indeed if c is large enough, the probability that a realization of the RW does not satisfy the SV property (1631 
decays as n~^. This statement can be easily shown, as in |30II33| . by using that for fixed j, Pr(|Xj| > x) < 
2exv(—x^/(2j) ) [th e so called Chernoff’s bound, see |49| p. 12[. Although the realizations of the RW that 


do not satisfy ( |163| l 
algorithm turns out 


necessitates a large number of probes ~ n, their contribution to the average cost of the 
to be negligible as they occur with a very small probability oc 1/n. On the other hand, as 
we shall see below, it is relatively easy to find the maximum of a RW which satisfies the “SV” property. 

The algorithm proposed by Odlyzko consists in two steps: 


• In a first stage, one searches a good estimate M*of M„. This is done by probing Xn, X 2 N, X^n,--- where 
N = [yVilognJ, where [xj denot es th e largest integer not larger than x. If the algorithm finds, here or later, 
a violation of the SV inequality (1631, one has to probe all the positions of the RW (b ut th is happens very 
rarely). We denote by M' = max Xq, Xn, X 2 N, X 3 N, ...} < Mn- If the RW satisfies SV (1631, then 


Indeed, if we denote by 


Mn — M' < cyjN log n = log 1 
such that Xk„,^^N < Mn < X(k 


(164) 

JV>^(;=max + l)Jv) < 

_ +i)iv)- As we 

discuss it below, it turns out that this estimate M' of Mn (1641 is however not precise enough for the forthcoming 
steps of the algorithm. It is indeed necessary to scan the neighborhood of the large X^jv’s on a finer window. 
If for some integer r one finds 


,,+i)jv then Mn — ma.x{Xk 
c^/IVlogrz, which follows from (1631, and which implies (1641 as A7 ^ max(7fic^^,^iv, 


XrN > M' — log 7 


(165) 


we probe XrN±jK, j = 1,2, ...I N/K\ ■ K ^ If the RW has SV, then any k with Xk = Mn must be as 

close as of a rN for some r for which (|165|) is true. We now denote M* the maximum of all probes found until 


now. Because we scan with intervals < log n around the maximum, the SV inequality (1631 give 


0 < Mn — M* < cy/EEogE < ' 


, 1/6 


(166) 
























33 


One can prove m that the average cost of this first phase of the algorithm is of order 0(^/n/ log n) negligible 
compared to the cost of the second phase, that we now describe, and which is of order 0{^/n). 


• With this estimate M* of the actual maximum M„, the second phase will eventually find in a number 
of probes that is of order 0{^/n), which is the leading contribution to the cost of this algorithm. To do this, 
we will scan the sample path from left to right as follows. We introduce m the index of the RW position Xm 
which is currently probed by the algorithm. We start with m = 0 and we denote by the greatest position 
probed so far by the algorithm including M*. At each step of this phase, two cases may occur: 


(i) 

(ii) 


If — Xm < the algorithm will probe the right neighbor of Xm and m is incremented by 1, 

m —>■ m + 1. 

If — Xm > this means that the algorithm is still far from the maximum, because we know that 

Mn — M* < n^^®. In this case, the immediate vicinity of Xm does not need to be explored and the strategy 
is to jump from Xm to Xm+k, where k is still to be determined. In order to be sure that the RW does not 


exceed M* between m and m + k, we must have in mind the upper envelope of the RW on that interval 
[m, m + k] (see Fig. 111. Hence we impose the following bound 


k < 2{M* - Xm) + {Xm - Xm+k) ■ 


(167) 


The first term in the right hand side of this inequality (1671, 2(AT^ — X m), is larger than , while the 
second term, {Xm — Xm+k) is boun ded by cy/klogn, thanks to SV (|l63|l - as stated above, if Xm+k — 


does not satisfy the SV inequality (1631, we abort this approach and probe every position. Hence we can 
choose k slightly smaller than 2{M* — Xm)- If m + k > n, we probe and stop. When the full path has 
been scanned, the maximum of the RW has been found by the algorithm. 


For a RW which satisfies SV ( |163| l, one can show |30| that the major contribution to the cost of the 
algorithm is when M* — Xm > Indeed, one can show that the contributions of the probes of the type 

(i) to the cost of the algorithm is of the order 0{n^^^). In fact, one can show that if the estimate M* of 
is such that M„ — M* < n“ then the cost of these contributions is of order 0{n^°‘). If we want that the cost 
of this part of the algorithm to be smaller than the cost of the last one, which is of order Oi^/n) . the n this 
requires 2a < 1/2, for instance 2a = 1/3, hence the choice a = 1/6 made by Odlyzko 1301 [see Eq. (|166|)|. 

The step size k is slightly smaller than 2(M„ — Xm) and we need only one probe to cont rol t he k positions 
between m and m + k. Since k can be interpreted as the velocity of the algorithm [see Eq. ( |159| l[, the average 
cost of the algorithm is, at leading order when n goes to inhnity, (C'(Od)) given by 


(Cl(Od)) 



1 

- Ai + 1 




(168) 


where we recall that (...) denotes an average over the different realizations of the RW XiS. When n goes to 
infinity, the RW becomes a BM and 


g(Od) ^ r dr 

n-s-oo 2 Jq a;max “ x{t) ’ 

as described in the text in (|^. 


(169) 


C.3 Odlyzko’s algorithm for the Bridge 


It is easy to check that the arguments presented above can be easily transposed to the case of a random walk 
bridge. In particular, given that the bridge is pinned at both extremities Ao,s = = 0, it s variations are 

typically smaller than the one of the free walk and hence the property of “slow variations” (1631, which plays a 


crucial role in this algorithm, would follow naturally. Therefore we conjecture that Odlyzko's algorithm w ould 
be the optimal one to Hnd the maximum Mn,B and its cost would be given by {l/2)T^-_i{t) given in Eq. (131. 


D Some useful integrals involving confluent hypergeometric functions relevant for the 
case V{x) — Ijx 

D.l An integral involving a single confluent hypergeometric function 

For the analysis of the functional Ta--i{t) [se e Eq . ( |127| )[, a useful integral involving the confluent hypergeo¬ 
metric function U{a,2,z) is the following (see |5U| as well as Mathematica): 













34 


PC 

(p{s) = 2^^^xr{x/V^) / 

^0 


— \/^y 


yU [1+ -^,2,2V^y^ dy 


= (— 271 A CSC ( —7=^ — 


V2s 


\V^ 


XHi-^-^^ 
2 2V^ 


)+ah(- 


2V^ 


(170) 

(171) 


where csc(x) = 1/sina; and H{x) are harmonic numbers, H{x) = 'ip{x) + 7 b where 'xp{x) = F'(x)/r(x) is the 
di-gamma function and 7 b the Euler constant. The function H{x) admits the following series expansion 


H{x) = Y.^-lY^{j + 2)x^+^ 

7 = 0 


(172) 


where ^(a:) is the Riemann ze ta fu nction. By combining (1701, together with (172 ) one arrives straightforwardly 
at the formula given in Eq. (11281) in the text. 


D.2 An integral involving the product of two confluent hypergeometric functions 


To compute the ampl itudes cb such that the functions (f)E{x) in ( |138[ ) with ds = 0 satis fy the orthogonality 
condition in Eq. (139) we used the following relation, derived by Landau and Lifshitz m (see formula (f.9) m 
Appendix f): 


POO 

J = / i-Fi(a, 7, kz)iFi{a ^ 7^, k'z)dz 

Jo 

= r(7)A“+“'-^(A - fc)-“(A - fc')"“' 2A1 (^a, a', 7, 


kk' 

{X-k){X-k')J ’ 


(173) 


where 2 ^ 1 ( 0 , q', 7 , 2 ;) is a generalized hypergeometric series. Such integrals (1731 arise naturall y in t he study 
of certain matrix elements of quantum Hamiltonian involving Coulomb interactions. In our case (138), one has 
= 1 — is/x/A, a' = 1 — isjy/W, 7 = 7 ' = 2, fc = 2iV2E. k' = 2iV2E' and A = 2\/2s. Hence the desired 


formula In our case (139) can b e obt ained by differentiating (1731 once wrt A and analyzing in detail the limit 
fc —>■ fc' o f the resulting formula ( |173[ |. These somewhat cumbersome manipulations yield the expression for cb 
given in (1401. 


E Numerical simulations of constrained Brownian motion 

In this appendix, we describe the algorith ms th at we have used here to simulate various constrained Brownian 
motions. We refer the interested reader to [52| for an extended discussion of these algorithms. 


E.l Brownian motion 

In order to simulate a Brownian motion x{t), we consider the discrete random walk 


fXo = 0 

\w = + [l,A] 


(174) 


where r/i’s are identical and independent Gaussian standard variables of variance unity. When N goes to 
infinity, A^iv] x{t) where x{t) is a Brownian motion, with r G [0,1] : x(t) = C(r), where ^(r) is a Gaussian 
white noise (C(x)C(t')} = S{t — t')- This is the building block (174), to simulate different constrained Brownian 
motions. 


//generation of Brownian Motion 

void BM ( int N, double *X, gsl_rng * r) 

{ 

//r is the ’seeds’ of the random number generator, 
int i; 

X[0]=0; 

for(i=l; i<N; i++) 

f 

X[i]=X[i-1]+gsl_ran_gaussian (r, l)/sqrt((double)N); 

} 

} 
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(a) Brownian motion 
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(c) Brownian excursion 
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(e) Reflected Brownian motion 
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(b) Brownian bridge 
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(d) Brownian meander 
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0 1 

(f) Reflected Brownian bridge 





Fig. 13 Example of different constrained Brownian motions studied in the present paper. 


E.2 Brownian bridge 

For a Brownian bridge xb{t), which is a Brownian motion starting and ending at the origin 3:s(0) = a:s(l) = 0, 
we use the identity x{t) — tx{1) = Xb{x) 


R -Xiv, [0,iV], 


(175) 


where XiS are generated by (174l. One can show that Rtjv] converges to a Brownian bridge xb{t). 


//generation of Brownian Bridge 

void BB ( int N, double *X, gsl_rng * r) 

{ 

BM(N,X,r); 
int i; 

for(i=l; i<N; i++) 

f 

X[i]=X[i] -(double)i/(N-l)*X[N-l] ; 

> 

} 


E.3 Brownian excursion 


For a Brownian excursion xe{t), which is a Brownian motion that starts and ends at the origin xb(0) = xe{1) = 
0 and staying positive in the interval [ 0 , 1 ], we use the identity \/[xb,i{t)\'^ + [a;B,2('r)]^ + [* 3 , 3 ( 1 ")]^ = xe{t) 
where xb,i,xb,2 and xb,3 are three independent Brownian bridges [SSlEll. Hence we simulate 


E. = i e [o,iv] 


(176) 


where yi,i,Y 2 ,i and are three independent realisations of (1751. FfrAr] converges to a Brownian excursion 
Xe{t). 

//generation of Brownian Excursion 
void BE ( int N, double *X, gsl_rng * r) 

{ 

double XI[N]; 
double X2[N]; 
double X3[N]; 
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BB(N,Xl,r) ; 

BB(N.X2,r); 

BB(N,X3,r); 
int i; 

for(i=l; i<N; i++) 

X [i] =sqrt (XI [i] *X1 [i] +X2 [i] *X2 [i] +X3 [i] *X3 [i] ) ; 

} 


E.4 Brownian meander 


For a Brownian meander XMeir), a Brownian motion which begins at the origin and stays positive on [0,1], 
one can show that the PDF of its Hnal position xf > 0 at time 1 is p{xf) = xfsT^^^"^ ■ One can then use the 
following representation of the meander ending at xf'- [53ll54| \/[xb,i{t)Y + [xB,2iT)Y + [xB,3{r) + txfY = 
xm^{t) where xb,i,xb.2 and xb.3 are three independent Brownian bridges and xf is a random variable drawn 
from p[xf) = XFeT^^^"^ ■ Hence the Brownian meander XMe('r) can be generated numerically as 




N 


i £ [0, N] 


(177) 


where Yi, Y 2 and Y 3 are three independent realizations of (1751, where / > 0 is a random variable, whose PDF 
is given by p{f) = /e“^ converges to a Brownian meander XMe('r). 


//generation of Brownian Meander 
void BMe ( int N, double *X, gsl_rng * r) 
{ 


double f=gsl_ran_rayleigh (r,l); 

double XI [N]; 
double X2 [N]; 
double X3[N]; 

BB(N,Xl,r); 

BB(N.X2,r); 

BB(N.X3,r); 
int i; 

for(i=0; i<N; i++) 

X [i] =sqrt ( (XI [i] +f* (double) i/ (N-1) ) * (XI [i] +f* (double) i/ (N-1) )+X2 [i] *X2 [i] +X3 [i] *X3 [i] ) ; 

> 

} 
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